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Abstract 

Considering the large value of the top Yukawa coupling, we investigate the lead- 
ing one- loop Yukawa electroweak corrections that can be induced by the top quark in 
a process such as Higgs production in association with a tagged bottom anti-bottom 
pair at the LHC. At NLO these contributions are found to be small at the LHC 
both for the total cross section and for the distributions. In the limit of vanishing 
bottom Yukawa coupling where the LO contribution vanishes, the process can still 
be induced at one-loop through the top quark transition. Though this contribution 
which can be counted as part of the NNLO correction is small for Higgs masses 
around 120GeV, it quickly picks up for higher Higgs masses. This contribution rep- 
resents the rescattering of the top quarks and their decay into VK's leading to Higgs 
production through WW fusion. 



1 Introduction 



The most important goal of the LHC is the discovery of the Higgs and the concomi- 
tant study of the mechanism of electroweak symmetry breaking. Especially if no new 
phenomenon is unravelled through the direct production of new particles, the study of 
the Higgs properties such as its self-couplings and couplings to the other particles of the 
standard model (SM) will be crucial in order to establish the nature of the scalar compo- 
nent of the model. In this respect most prominent couplings, in the SM, are the Higgs, 
the top, and to a much lesser degree the bottom, Yukawa couplings. The top Yukawa 
coupling is after all of the order of the strong QCD coupling and plays a crucial role in 
a variety of Higgs related issues. The dominant mechanism for Higgs production at the 
Large Hadron Collider (LHC) is the gluon fusion process, which incidentally is initiated 
through a top loop. Electroweak gauge boson fusion and W/ZH associated production 
P are also of importance. Higgs production associated with heavy quarks like the top 
or bottom quark is not considered as a discovery channel because of its small total cross 
section, the top suffering further from a complicated final state topology. However, if 
one wants to determine the bottom-Higgs Yukawa coupling, XbbH, then Higgs production 
associated with a bottom anti-bottom pair could provide a direct measurement of this 
coupling. In the minimal supersymmetric standard model, MSSM, the bottom Yukawa 
coupling is enhanced by a factor tan/5, the ratio of the vacuum expectation values of 
the two Higgs doublets. For high tan j3 and not too large Higgs masses this provides an 
important discovery channel for the supersymmetric Higgses. 

The next-to-leading order (NLO) QCD correction to pp —>■ bbH has been calculated by 
different groups relying on different formalisms. In a nut-shell, in the five-flavour scheme 
(5FNS)[2], [3], use is made of the bottom distribution function so that the process is ap- 
proximated (at leading order, LO) by the fusion bb — > H. This gives an approximation 
to the inclusive cross section dominated by the untagged low px outgoing b jets. If only 
one final b is tagged, the cross section is approximated by gb bH. The four flavour 
scheme (4FNS) has no b parton initiated process but is induced by gluon fusion gg bbH, 
with a very small contribution from the light quark initiated process qq Here 
again the largest contribution is due to low pt outgoing 6's which can be accounted for 
by gluon splitting into bb. The latter needs to be resummed and hence one recovers most 
of the 5FNS calculation while retaining the full kinematics of the reaction. QCD NLO 
corrections have been performed in both schemes [21 IH El E] and one has now reached a 
quite good agreement [7] . 

The 5FNS approach, which at leading order is a two-to-one process has allowed the com- 
putation of the NNLO QCD correction[8l [9] and very recently the electroweak/SUSY 
(supersymmetry) correction[10j to 66 — > 0, any of the neutral Higgs boson in the 
MSSM. SUSY QCD corrections have also been performed for gg bbh^Tl [12] where h 
is the lightest Higgs in the MSSM as well as to gb b6fi3\. 

In order to exploit this production mechanism to study the Higgs couplings to 6's, one must 

*In fact qq bbH is dominated by qq — > HZ* —^ bbH and does not vanish for vanishing bottom 
Yukawa couphng. However this contribution should be counted as ZH production and can be excluded 
by imposing an appropriate cut on the invariant mass of the bb pair. 
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identify the process and therefore one needs to tag both 6's, requiring somewhat large pt 
b. This reduces the cross section but gives much better signal over background ratio. For 
large px outgoing quarks one needs to rely on the 4FNS to properly reproduce the hight pt 
b quarks. The aim of this paper is to report on the calculation of the leading electroweak 
corrections to the exclusive bbH final state, meaning two 6's are detected. These leading 
electroweak corrections are triggered by top-charged Goldstone loops whereby, in effect, 
an external b quark turns into a top. This transition has a specific chiral structure whose 
dominant part is given by the top mass or, in terms of couplings, to the top Yukawa 
coupling. Considering that the latter is of the order of the QCD coupling constant, the 
corrections might be large. In fact, as we shall see, such type of transitions can trigger 
gg — > bbH even with vanishing XbbH in which case the process is generated solely at 
one-loop. We will quantify the effect of such contributions. 

This calculation belonging to the class of the 2 — >^ 3 processes at the LHC, we will 
also cover some technical issues pertaining to such calculations, like among other things 
the helicity amplitude method we use and the occurrence of vanishingly small (inverse) 
Gram determinant. This determinant occurs when solving the system of (linearly inde- 
pendent) tensor integrals in terms of the basis of scalar integrals. The Gram matrix is 
constructed out of the scalar products of the {N — 1) linearly independent momenta for a 
process with N external legs. The Gram determinant can vanish if the momenta of the set 
are, for example, for some exceptional point in phase space no longer linearly independent. 

In this paper we restrict ourselves to a Higgs mass in the range preferred by the latest 
electroweak data[14], in particular we confine the present study to Mh < ISOGeV. An- 
other reason for this choice is that, as we will briefly point out, as the Higgs mass increases 
the loop induced cross section increases and the loop integral starts showing instabilities. 
This we have identified as a Landau singularity which is a pinch singularity of the loop 
integral. This has an interesting physical origin: the rescattering of on-shell top quarks 
into W bosons, giving rise to W boson fusion into Higgs. We leave this important issue 
to another study though. 

The plan of the paper is as follows. In the next section we present some general con- 
siderations concerning the properties and structure of the calculation we have performed. 
We first briefly review the tree-level LO amplitude and highlight some symmetries of the 
helicity amplitudes. These symmetries are maintained by QCD corrections but not by 
the electroweak corrections we are studying. We then discuss the leading approximation 
as given by the insertions of the top-bottom-Goldstone Yukawa vertex. We classify the 
contributions into three classes in the cases of the NLO correction as well as the contribu- 
tions that survive at one-loop even for XbbH = 0. In section 3 we give our renormalisation 
scheme and discuss the inclusion of a top/Higgs Yukawa enhanced contribution which can 
be considered as a universal correction to Higgs processes related to the Higgs wave func- 
tion renormalisation and the renormalisation of the vacuum expectation value. Section 4 
gives an overview of some calculational details in particular how the calculation is organ- 
ised. Discussion on the loop integrals, the appearance of spurious instabilities related to 
vanishing Gram determinants and how these are cured depending on how the phase space 
integration is carried out is also presented. We also discuss in this section how we checked 
our results through ultraviolet finiteness and gauge invariance. Section 5 presents and 
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discusses the numerical results we find for the total cross section and various distributions 
both at the level of the NLO electroweak correction as well as the one-loop contribution 
that survives in the limit of vanishing XbbH- Details about the helicity amplitude method 
we used as well as the optimisation of the code are presented in two appendices. 

2 General considerations 

Before discussing the details of the calculation it is educative to expose some key features 
that appear when one considers the electroweak corrections at one-loop compared to 
the structure we have at tree-level or even the structure that emerges from QCD loop 
calculations. In particular the helicity structure is quite telling. So let us set our definition 
first. The process we consider is g{pi, Xi)+g{p2, X2) — ^ ^(ps, A3) + 6(p4, \4)+H{p^). Aj = ± 
with i = 1,2,3,4 are the helicities of the gluons, the bottom and anti-bottom while pi 
are the momenta of particles. The corresponding helicity amplitude will be denoted by 

2.1 Leading order considerations 



g(Pr ^) b(P3, A3) 




Figure 1: All the eight Feynman diagrams can be obtained by inserting the Higgs line to 
all possible positions in the bottom line. 

At tree- level, see Fig. [T]for the contributing diagrams, the Higgs can only attach to the 
6-quark and therefore each diagram, and hence the total amplitude, is proportional to the 
Higgs coupling to bb, XbbH- Compared to the gluon coupling this scalar coupling breaks 
chirality. These features remain unchanged when we consider QCD corrections. Moreover 
the QCD coupling and the Higgs coupling are parity conserving which allows to relate 
the state with helicities (Ai, A2; A3, A4) to the one with (— Ai, — A2; — A3, — A4) therefore 
cutting by half the number of helicity amplitudes to calculate. With our conventions for 
the definition of the helicity states, see Appendix [XJ parity conservation for the tree-level 
helicity amplitude gives 

-^0(^-^1) —-^2; "-^3) "-^4) = -^3-^4-^0 (-^1) -^2; -^3) -^4)*- (1) 

This can be generalised at higher order in QCD with due care of possible absorptive parts 
in taking complex conjugation. 

The number of contributing helicity amplitudes can be reduced even further at the 
leading order, in fact halved again, in the limit where one neglects the mass of the 6-quark 
that originates from the 6-quark spinors and therefore from the b quark propagators. We 
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should in this case consider the XbbH as an independent couphng, intimately related to the 
model of symmetry breaking. In this case chirality and helicity arguments are the same, 
the b and b must have opposite helicities for the leading order amplitudes and hence only 
^o(Ai, A2; A, —A) remain non zero. In this limit, this means that only a string containing 
an even number of Dirac 7 matrices, which we will label in general as r*^™"^ as opposed to 
-podd ^ string with an odd number of 7's, can contribute. 

In the general case and reinstating the b mass, we may write the helicity amplitudes as 

^(Ai,A2;A3,A4) = u{Xs){rz:i + n'X)<^^) 

= 5x„-x, (^^^^'^ + rribA'""') + 6x,m {■^°'''' + ^b^^™'^) . (2) 

The label '^^^^ in A'^^'^^ and are the contributions of r^^*^"^ to the amplitude and 

likewise for This way of writing shows that mb originates from the mass insertion 
coming from the massive spinors and are responsible for chirality flip. In the limit — 0, 
-'^aTa2 ■'^Ai'a2 contribute to different independent helicity amplitudes. In general F''™'^ 
and r°*^'^ differ by a (fermion) mass insertion. In fact r°'^'^ is proportional to a fermion 
mass insertion from a propagator. At leading order the mass insertion is naturally rrib, 
such that r°'^'^ is 0{mb)- This shows that at leading order, corrections from mi, = to 
the total cross section are of order 0{ml). Of course there might be some enhancement 
of the 0{mf) terms if one remembers that the cross section can bring about terms of 
order m1/{p\)^. However, in our calculation where we require the 6's to be observed 
hence requiring a cut, the effect will be minimal. With = 4.62GeV, the effect 
of neglecting nih is that the cross section is increased by 3.7% for Ip^^''] > 20GeV and 
1.1% for Ip^"! > 50GeV. At one-loop, the chiral structure of the weak interaction and the 
contribution of the top change many of the characteristics that we have just discussed for 
the tree-level. 

2.2 New electroweak Yukawa-type contributions, novel charac- 
teristics 




Figure 2: Sample of one-loop diagrams related to the Yukawa interaction in the SM. xw 
represents the charged Goldstone boson. 



Indeed, look at the two contributions arising from the one loop electroweak corrections 
given in Fig. [21 Now the Higgs can attach to the top or to the W. Therefore these 
contributions do not vanish in the limit XbbH = 0. Because now the fermion loop is a top 
loop, the mass insertion in what we called F"*^*^ is proportional to the top mass and is not 
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negligible. In fact the diagrams in Fig. [2] show the charged Goldstone boson in the loop. 
The latter triggers a t — > bxw transition whose dominant coupling is proportional to the 
Yukawa coupling of the top. We will in fact be working in the approximation of keeping 
only the Yukawa couplings. This reduces the number of diagrams and if working in the 
Feynman gauge as we do in this computation, only the Goldstone contributions survive. 
The neutral Goldstone bosons can only contribute corrections of order A^. We will neglect 
these 0{Xl) contributions at the amplitude level. However the order 0{Xb) corrections 
will be kept. All the corrections are then triggered by t — *• bxw and apart from the QCD 
g ^ bb vertex, only the Yukawa vertices shown in Fig. [3] below are needed to build up the 
full set of electroweak corrections. 




Figure 3: Relevant vertices appearing at one loop, eu = Xb/Xt, v is the vacuum expectation 
value and X is the Higgs self- coupling, related to the Higgs mass in the Standard Model. 
Pl,r = (It 

Note that in models outside the Standard Model, the Higgs coupling to the fermion 
/, XffH, can involve other parameters beside the corresponding Yukawa coupling A/. 
The Higgs coupling to the charged Goldstone involves the Higgs self-coupling or Yukawa 
coupling of the Higgs, A = Mfj/2v'^ proportional to the square of the Higgs mass. The 
latter can be large for large Higgs masses. These considerations allow to classify the 
contributions into three gauge invariant classes. 

2.3 Three classes of diagrams and the chiral structure at one- 
loop 

All the one-loop diagrams are classified into three gauge invariant groups as displayed 
in Fig. m The Higgs couples to the bottom quark in the first group (Fig. H^), to the 
top quark in the second group (Fig. |3)d) and to the charged Goldstone boson in the 
third group (Fig. Ht). As shown in Fig. H] each class can be efficiently reconstructed 
from the one-loop vertex bbH, depending on which leg one attaches the Higgs, by then 
grafting the gluons in all possible ways. We have also checked explicitly that each class 
with its counterterms, see below, constitutes a QCD gauge invariant subset. Note that 
these three contributions depend on different combinations of independent couplings and 
therefore constitute independent sets. 

The chiral structure t — > bxw impacts directly on the structure of the helicity ampli- 
tudes at one-loop. The split of each contribution according to F''™'^ and F"*^*^, see Eq. [2] 
will turn out to be useful and will indicate which helicity amplitude can be enhanced by 
which Yukawa coupling at one-loop. We show only one example in class (6) of Fig. |H It 
is straight forward to carry the same analysis for all other diagrams. We choose the first 



5 



Xw 



(b) 



V _ X 

Xw 



■ H 



\H(P5 






»- 




t , 


•^,Xw 




g(P2) - 


Xw 


b(P4) 



■000000 



'■0W0"0" 



*5 

\ 


I 


\ 

\ 

\ 

\ 


\ 




> N 

\ 




"1 ^ — 








— »■ 









''0000000' 



^0000000 



'"OOCCCOO" 



•000000 



'■ffOCDWO" 



' OOOOO W f — >• 



'■ffOOOOOO" 



■"0000000 



Figure 4: yl// i/ie diagrams in each groups can be obtained by inserting the two gluon lines 
or one triple gluon vertex (not shown) to all possible positions in the generic bottom line, 
which is the first diagram on the left. We have checked the number of diagrams through 
Grace-looplH 



diagram in group (b) in Fig. |H For clarity we will here take mb = 0, we have already 
shown how m;, insertions are taken into account, see Eq. [21 Leaving aside the colour part 
which can always be factorised out (see Appendix [B]) and the strong coupling constant, 
we write explicitly the contribution of this diagram as 



Ai(Ai, Aa; A3, A4) = \ttH>Hu{^-i.^P-i)i{K,Pi)^Ci,i^^{\2,P2)u{\i,Pi). (3) 

Cbi is the Yukawa vertex correction. In D-dimension, with q the integration variable, the 
momenta as defined in Fig. [1] with pij = Pi + Pj and Pij = Pj — Pi we have 

r = f {Pr - £btPL){mt + 4 + pi3){mt + 4 - MjPL - BbtPR) ... 

y (2vr)^z (M^ - q2)[m^ - (q + p,,nm^ - (q - p,,y] ' ^> 

The numerator of the integrand of (HI), neglecting terms of 0{Xl), can be re-arranged 
such as 

Ai(Ai, A2; A3, A4) -Ebt [ml + - M) +mPR (2^+^13 - 1^24) • (5) 

^ . ' ^ V ' 

pcvcn podd 

This shows explicitly that F"'^'^ structures with a specific chirality, Pr, can indeed be 
generated. They do not vanish as Xbth 0. The even one-loop structures on the other 



6 



hand are 0{\h). The structure in class (c), Higgs radiation off the charged Goldstones, 
is the same. For class (a), radiation off the 6-quark, the structure of the correction is 
different, the odd part is suppressed and receives an 0{Xb) correction. To summarise, 
with rrib = 0, making explicit the Yukawa couplings and the chiral structure if any, for 
example Pr, that characterise each class and comparing to the leading order, one has 





peven 


podd 


tree-level 
(a) 
(b) 
(c) 


Xb^t^ttH 
^b^t^XXH 




^tKxH, (Pr) 



We have kept X/jh and A/ separate to show how the structures may change in the 
MSSM for example and also why just by inspecting the couplings we can differentiate 
between the three classes. We clearly see that all one-loop F'^™" contributions vanish 
in the limit Af, = or Xf,bH = 0. On the other hand this is not the case for the one- 
loop F°*^'^ contribution belonging to class (6) and (c). However for these contributions 
to interfere with the tree-level LO contribution requires a chirality flip through a rrib 
insertion. Therefore in the SM for example, the NLO cross section is necessarily of order 
ml, like the LO, with corrections proportional to the top Yukawa coupling for example. 
On the other hand, in the limit of XbbH = 0, the tree level vanishes but gg bbH still goes 
with an amplitude of order g'^X'^XuH or g^X^X^^n- For XbbH 7^ these contributions should 
be considered as part of the NNLO "corrections" however they do not vanish in the limit 
rub — *■ (or XbbH = 0) while the tree level does. These contributions can be important 
and we will therefore study their effects. For these contributions at the "NNLO" we can 
set nib = 0. 

The classification in terms of structures as we have done makes clear also that the novel 
one-loop induced F°'^'^ contributions must be ultraviolet finite. This is not necessarily the 
case of the structures where counterterms to the tree-level structures are needed 

through renormalisation to which we now turn. 

3 Renormalisation 

We use an on-shell (OS) renormalisation scheme exactly along the lines described in|15]. 
Ultraviolet divergences are regularised through dimensional regularisation. In our approx- 
imation we only need to renormalise the vertices bbg and bbH as well as the bottom mass, 
rUb- For the bbg vertex, from the point of view of the corrections we are carrying, only 

{-) ^ , (-) 

wave function renormalisation for the b/b field is required: bL,R — ^ (1 + ^^bLR) bL,R- 
5Zl^^^ can be taken real, see[l5]. The counterterm to mj,, 5mb and the wave function 
renormalisation for the b/b are set by imposing the usual conditions for pole position and 
residue on the renormalised bottom propagator. In terms of the self-energy correction 
Sbb(g) with momentum g[T5]: 

nq^) = Kr + K^4 + K^^i-f^. (6) 
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This translates into 

Snih = Re(^mbK^{ml) + Ki{ml)^, 
bZl'^ = ^Re(^K5^{ml) - K^{ml)^ - mh-^(mbReK^{q'^) + ReKi{q^)^\g2^^2, 

^^bR = -^Re(^K^^{ml) + K^iml)^ - mh-^(mhReK^{q'^) + ReKi{q^)^ |g2=„2.(7) 

We calculate the coefficients Ki^^^^^ of the bottom self-energy in the same spirit we cal- 
culate the other one-loop corrections, i.e. only through the t bxw transition, see the 
first diagram of class (a) in Fig. HI We get. 



Ki{q') = (Cuv-Fo{mt,Mw,q')), 
K,{q^) = -K,,{q^) = ^ [Cuv - 2Fi(m„ M^., q^)) 
with C^y = i - 7£; + ln47r ,£) = 4 - 2e, 

The reason we get K^{q^) = —K^^i^q"^) is due to the particular chiral structure of the 
t bxw loop insertion. In particular for = 0, one recovers that these corrections only 
contribute to 6Z}^'^ and not dZ^'^. 

Oh bji 

The counterterms needed to renormalise the bbH vertex are 6mh, dzV"^ as well as the 

1/2 

Higgs wave function renormalisation 5Z^ and the counterterm to the vacuum expectation 
value, v,bv. Indeed we have = XbbH{^^ + ^Z^'^ +6Z^'^ +{5Z^ —5v)). The t — > bxw 

-1 /o 

loop insertion does not contribute to 5Z^ (which originates from the Higgs self energy 
two-point function) nor to the renormalisation of the vacuum expectation value. On 

1 /2 

the other hand {5Z^ — 5v) can be seen as a universal correction to Higgs production 
processes. We will include this correction as it has potentially large contributions scaling 
like and A which fall into the category of the corrections we are seeking. Within the 

1 /2 

calculation we have performed this means that the combination {5Z^ — 6v) must be 
finite. Indeed, we find 

1 r SA^ / 

6Z]I' = -^—Rel-^(Cuv-Fo{m^,m„Mfj) 

- M^Goimt, mt, Ml) + AmlG^{mu M|)) 

- ^ (9^0(^1^^, Mu, Ml) + 2GoiMw, Mw, M^) + GoiMz, Mz, M|,)) | , 
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Sv = -—Rel-^(Cuv-2F,{m,,mt,M^)J 

- x(^Fo{Mh,Mw,M^) - F,{Mh,Mw,M^) - UnMjj^ 

- T (^o(mt, mt, Ml) - 2Fi(m,, m„ M^)) + X{Fo{Mh, Mz, M|) 

- Fi(Mh, Mz, M|) - Fo{Mh, Mw, M^) + Fi{Mh. Mw, M^))) } , 

d — x^xfl — 

G„(mi, ms, g^) = g^— F„(mi, ma, g^) = gM c^a;- 2 71 TT' 

ag^ Jo ~ xjmf + — x{l — x)q^ 

(9) 

1 /2 

which shows that {SZ^ — 5v) is finite. 

In the actual calculation, the counter term 5^^^ belongs to class (a) in the classification 
of Fig m This makes class (a) finite. The counterterm we associate to class (6) is the 
part of 5hbH from the t — >• hxw loops and therefore does not include what we termed the 
universal Higgs correction, i.e does not include the contribution {SZ^^"^ — 5v). This is 
sufficient to make class (6) finite. In our approach (c) is finite without the addition of a 
counterterm. We will keep the {SzH"^ — 6v) contribution separate from the contributions 
in classes (a), (6), (c). We will of course include it in the final result. 

4 Calculation details 

We have written two independent codes. In the first one we set m^ = in all propagators 
and other spinors that emerge from the helicity formalism we follow. In this limit, the 
helicity formalism is very much simplified and the expression quite compact. This code 
is in fact subdivided in two separate sub-codes. One sub-code is generated for the "even" 
part (constituted by the F*^™" contributions, see Eq. [21) and the other by the "odd" part. 
We also generate a completely independent code for the case m^ 7^ where in particular 
we use the helicity formalism with massive fermions. Details of the helicity formalism 
that we use are given in Appendix [XI 

The steps that go into writing these codes are the following. In the first stage, we 
use Form[T6] to generate expressions for the tree level and one loop helicity amplitudes. 
Each helicity amplitude is written in terms of Lorentz invariants, scalar spinor functions 
{A, B,C)x-x defined in Appendix [XI and the Passarino- Veltman fU] tensor functions T^J 
for a tensor of rank M for A^-point function. We have also sought to write the contribution 
of each amplitude as a product of different structures or blocks that reappear for different 
graphs and contributions. For example colour factorisation is implemented, this further 
allows to rearrange the amplitude into an Abelian part and a non-Abelian part which will 
not interfere with each other at the matrix element squared level. The helicity information 
is contained in a set of basic blocks for further optimisation. Another set of blocks pertains 
to the loop integrals and other elements. The factorisation of the full amplitude in terms 
of independent building blocks is easily processed within Form. These building blocks can 
still consist of long algebraic expressions which can be efficiently abbreviated into compact 
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variables with the help of a Perl script which also allows to convert the output of Form 
into the Fortran code ready for a numerical evaluation. More details on the connection 
between Form and Fortran as well as the optimisation we implemented in the codes can 
be found in Appendix [Bl 



4.1 Loop integrals, Gram determinants and phase space inte- 
grals 

The highest rank M of the Passarino-Veltman tensor functions with M < N that we 
encounter in our calculation is M = 4 and is associated to a pentagon graph, = 5. We 
use the library LoopTools[18] to calculate all the tensorial one loop integrals as well as 
the scalar integrals, this means that we leave it completely to LoopTools to perform the 
reduction of the tensor integrals to the basis of the scalar integrals. In order to obtain the 
cross section one needs to perform the phase-space integration and convolution over the 
gluon distribution function (GDF), g{x,Q) with Q representing the factorisation scale. 
We have 

a{pp^bbH) = -i- / dxig{xi,Q) I dx2g{x2,Q) 
256 Jo Jo 

X i / ^^^\Aigg^bbH)mp,+p2-ps-p,-p,), 
F J 2e3 2e4 2e5 

(10) 

where 2i6~i^l^l^^ ^P^^ ^"^^ colour average factor and the flux factor is 1/F = 
l/(2nf2s^ with s = xiXas > {2mb + Mnf- 

The integration over the three body phase space and momentum fractions of the two 
initial gluons is done by using two "integrators": BASES |I9] and DADMUL[20j. BASES is a 
Monte Carlo that uses the importance sampling technique while DADMUL is based on the 
adaptive quadrature algorithm. The use of two different phase space integration routines 
helps control the accuracy of the results and helps detect possible instabilities. In fact 
some numerical instabilities in the phase space integration do occur when we use DADMUL 
but not when we use BASES which gives very stable results with small integration error, 
typically 0.08% for 10^ Monte Carlo points. For the range of Higgs masses we are studying 
in this paper, the instabilities that are detected with DADMUL were identified as spurious 
singularities having to do with vanishing Gram determinants for the three and four point 
tensorial functions calculated in LoopTools by using the Passarino-Veltman reduction 
methocll]. Because this problem always happens at the boundary of phase space, we can 
avoid it by imposing appropriate kinematic cuts in the final state. In our calculation, 
almost all zero Gram determinants disappear when we apply the cuts on the transverse 
momenta of the bottom quarks relevant for our situation, see section 15.11 for the choice 
of cuts. The remaining zero Gram determinants occur when the two bottom quarks or 
one bottom quark and the Higgs are produced in the same direction. Our solution, once 



tTlie reduction of the five point function using the method of Dcnner and Dittmaier [211 [22 which 
avoids the Gram determinant at this stage as implemented in LoopTools gives very stable results. 
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identified as spurious, was to discard these points by imposing some tiny cuts on the 
polar, 9, and relative azimuthal angles, of the outgoing 6-quarks, the value of the cuts 
is ^cut — I sin0^|cut = 10~^. DADMUL then produces the same result as BASES within the 
integration error. 

4.2 Checks on the results 

i) Ultraviolet finiteness: 

The final results must be ultraviolet (UV) finite. It means that they should be indepen- 
dent of the parameter Cuv defined in Eq. |H1 In our code this parameter is treated as a 
variable. The cancellation of Cuv has been carefully checked in our code. Upon varying 
the value of the parameter Cuv from Cuv = to Cuv = 10^, the results is stable within 
more than 9 digits using double precision. This check makes sure that the divergent part 
of the calculation is correct. The correctness of the finite part is also well checked in our 
code by confirming that each helicity configuration is QCD gauge invariant. 

ii) QCD gauge invariance: 

In the physical gauge we use, the QCD gauge invariance reflects the fact that the gluon is 
massless and has only two transverse polarisation components. In the helicity formalism 
that we use, the polarisation vector of the gluon of momentum p and helicity A is con- 
structed with the help of a reference vector q, see Appendix El for details. The polarisation 
vector is then labelled as e'^{p, X;q). A change of reference vector from q to q' amounts 
essentially to a gauge transformation (up to a phase) 

A; q') = d'^^^''''h^{p, A; q) + q)p>^ . (11) 

QCD gauge invariance in our case amounts to independence of the cross section in the 
choice of the reference vector, q. We have carefully checked that the numerical result for 
the norm of each helicity amplitude at various points in phase space is independent of the 
reference vectors say gi,2 for gluon 1 and 2, up to 12 digits using double precision. By 
default, our numerical evaluation is based on the use of qi^2 = (P2,Pi)- For the checks in 
the case of massive b quarks the result with the default choice qi^2 = (^2,^1) is compared 
with a random choice of qi^2, keeping away from vectors with excessively too small or too 
large components, see Appendix |X] for more details. 

iii) As stated earlier, the result based on the use of the massive quark helicity amplitude 
are checked against those with the independent code using the massless helicity amplitude 
by setting the mass of the b quark to zero. This is though just a consistency check. 

iv) At the level of integration over phase space and density functions we have used two 
integration routines and made sure that we obtain the same result once we have properly 
dealt with the spurious Gram determinant as we explained in section 14.11 

v) Moreover, our tree level results have been successfully checked against the results of 
CalcHEP^j. 
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5 Results 



5.1 Input parameters and kinematical cuts 

Our input parameters are a(0) = 1/137.03599911, Mw = 80.376 6GeV, Mz = 91.1876GeV, 
as{Mz) = 0.118, rub = 4.62GeV, mt = 174.0GeV with sw = ^/T^^M^JM^. The CKM 
parameter Vtb is set to be 1. We consider the case at the LHC where the center of mass 
energy of the two initial protons is ^/s = 14TeV. Neglecting the small light quark initiated 
contribution, we use CTEQ6L[21] for the GDF in the proton. The factorisation scale for 
the GDF and energy scale for the strong coupling constant are chosen to be Q = Mz for 
simplicity. 

As has been done in previous analyses [3, EB], for the exclusive bbH final state, we 
require the outgoing b and b to have high transverse momenta |p^^| > 20GeV and pseudo- 
rapidity \r]'''''\ < 2.5. These kinematical cuts reduce the total rate of the signal but also 
greatly reduce the QCD background. As pointed in [3] these cuts also stabilise the scale 
dependence of the QCD NLO corrections compared to the case where no cut is applied. 
In the following, these kinematical cuts are always applied unless otherwise stated. 

Talking of the NLO QCD scale uncertainty and before presenting our results, let us re- 
mind the reader of the size of the QCD corrections. Taking a renormalisation/factorisation 
scale as we take here at Mz, the QCD corrections in a scheme where the bottom Yukawa 
coupling is taken on-shell amount to ~ —22% for a Higgs mass of 120GeV. 



5.2 NLO EW correction with X^bH 7^ 




Figure 5: Left: the LO and NLO cross sections as functions of M^- Right: the rela- 
tive NLO EW correction normalized to tree level (Tlo- (c^)? (^)j (c) correspond to the 
three classes of diagrams as displayed in Fig. to which counterterms are added (see 
section\B^. [SZ]^^ — Sv) is the correction due to the universal correction contained in the 
renormalisation of the bbH vertex. "Total" refers to the total electroweak correction, of 
Yukawa type, at one-loop. 



The cross sections with two high-p^ bottom quarks at LO and NLO at the LHC are 
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displayed in Fig. [5] as a function of the Higgs mass. The NLO EW correction reduces the 
cross section by about 4% to 5% as the Higgs mass is varied from llOGeV to ISOGeV. 
The first conclusion to draw is that this correction is small if we compare it to the QCD 
correction or even to the QCD scale uncertainty. Considering that we have pointed to 
the fact that the contributions could be grouped into three gauge invariant classes that 
reflect the strengths of the Higgs coupling to the b, the t or its self-coupling, one can 
ask whether this is the result of some cancellation. It turns out not to be the case. All 
contributions are below 3%, see Fig. O Class (a) with a Higgs radiated from the bottom 
line is totally negligible ranging from —0.09% to —0.06%. We have failed in finding a 
good reason for the smallness of this contribution compared to the others. Those due 
to the Higgs self-coupling are below 1%. Radiation from the top contributes about —2% 
and is of the same order as the contribution of the universal correction. We had argued 
that the Yukawa corrections brought about by the top might be large. It seems that the 
mass of the top introduces also a large scale which can not be neglected compared to the 
effective energy of the hard process even for LHC energies. 




Figure 6: Effect of the NLO electroweak corrections on the pseudo-rapidity and transverse 
momentum distributions of the Higgs for = 120, ISOCeK The relative corrections 
dcTNLo/dcTLO — ^ is also shown. 



The NLO corrections are spread rather uniformly on all the distributions we have 
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looked at. We have chosen to show in Fig. [6] the effect on pseudo-rapidity and transverse 
momentum distributions of the Higgs for two cases = 120GeV and Mh = ISOGeV. As 
Fig. [H] shows the relative change in these two distributions is sensibly constant especially 
for Mh = 120GeV. For Mh = ISOGeV, the corrections are largest for p!f around 140GeV, 
however this is where the cross section is very small. A similar pattern, i.e. a constant 
change in the distributions, is observed for the bottom variables. 

5.3 EW correction in the limit of vanishing Xi,1)H 




MnlGeV] MufGeV] 



Figure 7: The one-loop induced cross section as a function of Mh in the limit of vanishing 
bottom-Higgs Yukawa coupling. The right panel shows the percentage contribution of this 
contribution relative to the tree level cross section calculated with X^hH 7^ 0. 

The cross section for Xh^n = can be induced at one-loop through the top loop. This 
"NNLO" contribution rises rather quickly as the Higgs mass increases even in the narrow 
range Mh = 110 — ISOGeV as can be seen in Fig. [71 Indeed relative to the tree level, the 
cross section with Mh = 120GeV amounts to 3% while for Mh = ISOGeV it has increased 
to as much as 17%. Going past Mh > 2M]y we encounter a Landau singularity [26] (a 
pinch singularity in the loop integral) from diagrams like the one depicted in Fig. [2] (right) 
with the Higgs being attached to the W^s or their Goldstone counterpart. It corresponds 
to a situation where all particles in the loop are resonating and can be interpreted as the 
production and decay of the tops into (longitudinal) W's with the later fusing to produce 
the Higgs. This leading Landau singularity is not integrable, at the level of the loop 
amplitude squared and must be regulated by the introduction of a width for the unstable 
particles. We leave this issue together with a general discussion of Landau singularities 
in such situations to another publication. 

Fig. [H] shows the pseudo-rapidity and transverse momentum distributions of the Higgs 
as well as the the px of the bottom for two cases Mh = 120GeV and Mh = ISOGeV 
in the limit of vanishing bottom-Higgs Yukawa coupling. These distributions are signif- 
icantly different from the ones we observed at tree-level (and with the electroweak NLO 
corrections), see Fig. [6l The Higgs prefers being produced at high value of transverse 
momentum, about ISOGeV. In the case of a Higgs with Mh = ISOGeV this contribution 
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Figure 8: The pseudo-rapidity of the Higgs and transverse momentum distributions of the 
Higgs and the bottom for Mh = 120, 150GeV arising from the purely one-loop contribution 
in the limit of vanishing LO (XbbH = 0/ Its relative percentage contribution da{\hbH = 
0)/daLo is also shown 
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can significantly distort the shape of the pff distribution for hight pif- with a " correction" 
of more than 70% over a rather large range. The distribution in the pt of the bottom is 
also very telling. The new contributions do not produce the bottom preferentially with 
low as the case of the LO contribution. 

6 Conclusions 

We have calculated the EW radiative corrections triggered by the Yukawa coupling of 
the top to the process pp —>■ bbH at the LHC through gluon fusion in the Standard 
Model. This process is triggered through Higgs radiation of the bottom quark with a small 
coupling proportional to the mass of the bottom. Yet in order to analyse this coupling, 
precision calculations that include both the QCD and electroweak corrections are needed. 
In this perspective, to identify the process one needs to tag both 6-jets. Our calculation is 
therefore conducted in this kinematical configuration. Inserting a top quark loop with a 
Yukawa transition of the type t bxw, Xw is the charged Goldstone, allows now the Higgs 
to be radiated from the top or from the Goldstone boson. The latter coupling represents 
the Higgs self-coupling and increases with the Higgs mass. The former, the top Yukawa 
coupling, is also large. As a consequence, the one-loop amplitude gg — » bbH no longer 
vanishes as the Higgs coupling to 6's does, like what occurs at leading order. We find that 
in the limit of vanishing XbbH, the one-loop induced electroweak process should be taken 
into account for Higgs masses larger than 140GeV or so. Indeed, though this contribution 
is quite modest for a Higgs mass of llOGeV it increases quite rapidly as the Higgs mass 
increases, reaching about 17% of the leading order value, calculated with = 4.62GeV, 
for Mh = ISOGeV. For these new corrections to interfere with the leading order requires 
helicity flip. Therefore at next-to-leading order in the Yukawa electroweak corrections, all 
corrections involve either a bottom mass insertion or a bottom Yukawa coupling. At the 
end the total Yukawa electroweak NLO contribution brings in a correction which is within 
the range —4% to —5% for Higgs masses in the range llOGeV < Mh < ISOGeV. They 
are therefore negligible compared to the NLO QCD correction and even the remaining 
QCD scale uncertainty. This modest effect translates also as an uniform rescaling of the 
distributions in the most interesting kinematical variables we have looked at (pseudo- 
rapidities and Pt of both 6-quarks and the Higgs). This is not the case of the one- loop 
induced contributions which survive in the limit of m;, ^ (and Xf,bH 0). Here the 
distributions for the Higgs masses where the corrections for the total cross section is 
large are drastically different from the LO distributions. A summary for the corrections 
including the NLO with Xbtn 7^ and the part of the NNLO counted as loop induced in 
the limit XbbH ^ is shown in Fig. [H 

The analysis we have performed in this paper does not cover Higgs masses over 150GeV 
and rests within the range of Higgs masses preferred by indirect precision measurements. 
In fact as the threshold for H —>■ WW opens up, important phenomena take place. Fore- 
most a Landau singularity, or a pinch singularity in some loop integrals, develops. This 
corresponds the rescattering of on-shell top quarks that decay to on-shell W with Higgs 
production via WW fusion. We leave this important issue to a forthcoming publication 
especially that the identification and handling of such singularities can be applied to other 
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Figure 9: Sew = Snlo + — - as a function of Mh- 



processes. In our case the singularity can be tamed by introducing the width of the un- 
stable particles. At NLO, for Mh = 2Mw for example, the wave function renormalisation 
of the Higgs, which involves the derivative of the two-point function Higgs self-energy, 
diverges. This can also be regulated by including the width of the W, see for example [27]. 

There is another contribution which does not vanish for vanishing XhtH and which 
contributes to gg — *■ bbH through a closed top quark loop. This contribution represents 
gg Hg* — > Hhh. We have not included this contribution in the present paper as we 
do not consider it to be a genuine bbH final state. This correction can be counted as 
belonging to the inclusive gg H process. The same line of reasoning has been argued 
in [7] . Nonetheless from the experimental point of view it would be interesting to include 
all these effects together with the NLO QCD corrections and the electroweak corrections 
that we have studied here. 
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Appendices 



A The helicity amplitude method 
A.l The method 

We use a combination of helicity amplitude methods as described in [281 129] to calculate 
the total cross section. In the following we only want to highlight some key features that 
were most useful for our calculation, for details of the method we refer to[2Hl[2n]- For our 
process g{pi, Xi) +g{p2, A2) — > ^(ps, A3) +b{p4, A4) +if(p5) where the particles are denoted 
by their momentum pi and helicity Aj we write the corresponding helicity amplitude as 
>A(Ai, A2; A3, A4). 



^(Ai,A2; A3,A4) = ef,{pi,\i;qi)eu{p2,X2;q2)M^''{\3,X4), 

-M^'^(A3,A4) = u{ps,X3)n{p,,X,). (A.l) 

r^'^ is a string of Dirac 7 matrices. These 7 matrices represent either interaction vertices 
or momenta from the fermion propagators. In our case the interaction vertices are the 
vectorial gluon vertices in which case they represent , the scalar Higgs vertex and at one- 
loop the pseudo-scalar Goldstone coupling. For the momenta, in our implementation we 
re-express them in terms of the independent external momenta Pi,P2,P3,P4- This applies 
also to the loop momenta after the reduction formalism of the tensor integrals has been 
performed. The first step in the idea of the helicity formalism we follow is to turn each of 
these 7 matrices (apart from the pseudo-scalar and the trivial scalar) into a combination 
of spinor function uu. We therefore transform our helicity amplitude into products of 
spinors such as the helicity amplitude could be written like a product u uu ...uu v with 
the possible insertion of 75's in the string. The different u, u, v in the string we have 
written have of course, in general, different arguments. Nonetheless one can turn each 
spinor product of two adjacent uu, etc into a complex number written in terms of the 
momenta in our problem as we will see. 
In the first step, for the momentum with pf = mf we use 

= u{pi, -)u{pi, -) + u{pi, +)u{pi, +) - rui. (A. 2) 

The polarization vector of the initial gluon i, e^{pi, A^; g^), is also first expressed in terms 
of spinors such as 



f . ^ u(pi, Ai)7^M(gi, Ai) 
e.(p.A,g.) = [4(^^.^^)]V2 ^ (A-3) 



where Qi is an arbitrary reference vector satisfying the following conditions 

g-=0, P^■q^^Q (A.4) 

Gauge invariance (transversality condition) requires that the cross sections are indepen- 
dent of the choice of the reference vector as we will see later. This acts as an important 
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check of the calculation, see later. It is not difficult to prove that the choice (1A.3P satisfies 
all the conditions for a transverse polarization vector. In particular, 

Pi.e{pi,\i) = 0, e{pi,\i).e{pi,Xi) = 0, 

€^{Pi,-\i) = e^{pi,\iy, e{pi,\i).e{pi,-\i) = -1, (A. 5) 

where the reference vector is not written down explicitly, i = 1,2 and no sum over i must 
be understood. Then for £^7^^ one uses the so-called Chisholm identity 

u{p, A)7^u(g, A)7^ = 2[u{p, -X)u{q, -A) + u{q, X)u{p, A)], (A.6) 

where all the spinors in Eq. IA.6I are for massless states in view of the lightlike condition 
on the reference frame vector and of course the momentum of the real gluon. 
With U{pi, Aj) representing either u{pi, Aj) or v{pi, Aj) one uses the general formulae 



U{pi, Xi)U{pj, Xj, 



Ax^XjiPhPj) + MiBx,x^{pi,pj) + MjCx,x,iPi,Pj] 
\/{Pi-ko){Pj-ko) 



V{Pi-ko){Pj-ko) 

where 

Mi = +mi if U{pi,Xi) = u{pi,Xi), 

Mi = -rrii if U{pi, Xi) = v{pi, Xi), 
Ax,x, = Sx,^x,Xi [{kQ.pi){ki.pj) - {ko.pj){ki.pi) - iXie^^p^k^^k'^p'^p'^) , 
Bx^x, = Sx,x,{ko.pj), Cx^x, =Sx,x,{ko.pi), (A.8) 

with /cq,! being auxiliary vectors such that k^ = 0, kf = —1 and fco-^i = 0. No sum over 
repeated indices must be understood. For instance, we can choose ko = (1,0,1,0) and 
ki = (0, 1,0,0). With this choice, it is obvious to see that the denominator in ( ]A.7ll can 
never vanish if the bottom mass is kept. If one would like to neglect rrib, that choice can 
bring p^.kQ or P4./C0 to zero in some cases. If this happens, one can tell the code to choose 
ko = (1, 0, —1, 0) instead of the above choice. In fact, that is what we did in our codes. 

In the case of spinors representing a massless state, the helicity formalism simplifies 
considerably. Only AxiXj is needed. Traditionally we introduce the C-numbers s{p, q) and 
t{p,q), 

s{p, q) = u{p, +)u{q, -) = A+_{p, q), t{p, q) = u{p, -)u{q, +) = -s{p, q)*. (A.9) 

These are the functions that appear in our code for the massless b quark. The massless 
case is also used when expressing the gluon polarisation vector to which we now turn. 



A. 2 Trans versality and gauge invar iance 

The reference vector used for the polarisation of the gluon can be changed at will. Chang- 
ing the reference vector from q to q' amounts essentially to a gauge transformation. Indeed 
one has [2H] 

A; q') = e^^(^' "^e'^lp. A; q) + g)p^ (A.IO) 
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where 



i(l>{q',q) 



P{q', q) 



nl/2 



q) t{p, q') 
t{p, q) s{p, q') 
2 t{q,q') 
[4(g'.p)]i/2 t{q,p)- 



(A.ll) 



Therefore up to the phase factor, the difference is contained in the momentum vector of 
the gluon. QCD gauge invariance for our process leads to the important identity 



|^(Ai,A2;A3,A4;gi,g2)|' = |^(Ai, A2; A3, A4; gl, gs)!', (A.12) 

as long as 2 satisfy the condition ( 1A.4I) . We have carefully checked that the numerical 
result for the norm of each helicity amplitude at various point in phase space is indepen- 
dent of the reference vectors gi,2 up to 12 digits using double precision. By default, our 
numerical evaluation is based on the use of gi^2 = (P2;Pi)- For the checks in the case of 
massive h quarks the result with gi 2 = {p2,Pi) is compared with the one using any gi^2 
such as the conditions (1A.4I) are obeyed. In the case of massless b quarks it is simplest to 
take gi,2 = {p3,Pi>- 

This check is a an important check on many ingredients that enter the calculation: 
the Dirac spinors, the gluon polarization vectors, the propagators, the Lorentz indices, 
the loop integrals. It has been used extensively in our numerical calculation. 



B Optimisation 

Each hehcity amplitude ^(Ai, A2; A3A4) = ^(A), a C-number, is calculated numerically 
in the Fortran code. The price to pay is that the number of helicity amplitudes to be 
calculated can be large, 16 in our case for the electroweak loop part. Some optimisation is 
necessary. The categorisation of the full set of diagrams into three gauge invariant classes 
as shown in section [23] is a first step. We have sought to write each diagram as a compact 
product of blocks and structures containing different properties of the amplitude. We 
write the amplitude according to a colour ordering pattern that defines three channels. 
The ordering is in a one-to-one correspondance with the three channels or diagrams shown 
in Fig. [H The T-type is the direct channel, the ?7-type is the crossed one obtained from 
the T-type by interchanging the two gluons and the S'-type is the one involving the triple 
gluon vertex. The helicity amplitude for each diagram can thus be represented as 

^(A)^'^'^ = CME{a, b)xCcx FFE x SME{Xi), (B.l) 

where 

• CME(a, b) is the colour matrix element, a, b are the colour indices of the two initial 
gluonq^l. The colour products can be (T^T''), (T^T'^) or [T'^,T^] corresponding to the 
3 T, U, S channels respectively 

"''Other colour indices of the bottom quarks are omitted here for simplicity 
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• Cc contains all the common coefficients like the strong coupling constant gg or factors 
common to all diagrams and amplitudes such as the normalisation factor entering 
the representation of the polarisation vector of the gluon, see Eq. IA.3I 

• FFE, form factor element, contains all the denominators of propagators, loop func- 
tions as well as various scalar products of external momenta {pi,P2,P3,P4} i-e. all 
the scalar objects which do not depend on the helicity Aj 

• SME{X), standard matrix element, is a product of the scalar spinor functions A\-\., 
BxiXj and Cx^\^ defined in Appendix IK\ 

For each channel, say ^(A)"^, the most complicated and time-consuming part is the 
FFE. That is why we want to factorise it out and put it in a common block so that 
in order to calculate all the 16 helicity configurations of ^(A)-'" we just need to calculate 
FFE once. This is done at every point in phase space. This kind of factorisation can be 
easily carried out in Form. 

SME{X) is also complicated because the bottom quark is massive and 75 occurs in the 
"helicity strings". Thus we have to optimize this part as well. The way we do it for all 
the 3 groups is as follows. In Form, we have to find out all the generic expressions of 
SME(\). There are 12 of them at tree level and 68 at one-loop if we choose qi^2 = P2,i 
for the reference vectors. For instance, 

SMEi = [u{Xs,p3)v{X4,p4)] X [e^{Xi,pi,p2)P4] x [e^{X2,P2,Pi)pl], 

= BMEi{X3,X^) X BME2{Xi) x BMEs{X2), (B.2) 

can be expressed in terms of 3 basic matrix elements (BME). Each BME occurs several 
times when calculating all the SME{X). The number of BME is 31. Each BME is 
written in terms of scalar spinor functions Ax^x^, Bx^x^, ^KXj- the SME or BME can 
be found and abbreviated in Form. As an alternative, we can use Perl for such an opera- 
tion. The Form output is converted directly into a Fortran code for numerical evaluation. 
Needless to say, all the abbreviations of SME or BME must be put in common blocks. 

To get the final result, we have to sum over all the channels. The grouping can be 
re-arranged in terms of an Abelian part and a non-Abelian part according to 

^(A) = ^(Af + ^(A)^ + ^(A)^, 

= {T^T^}^(A)^^"' + [T^TV(A)^^'"^ (B.3) 

where 

^(A)^^^' = \{A(xf + A(xf), 
^^^^TVAfed ^ j^(xf j^]^{A{Xf - A{Xf) (B.4) 

corresponding to the Abelian and Non-Abelian parts respectively. The amplitude squared 
then contains no interference term between the Abelian and Non-Abelian parts: 

I ^(A) |2= -i^ ('^ I ^(A)^"^' |2 +12 I ^(A)^^^^' A (B.5) 
where 2i6~i^l^l^^ spin- and colour- averaging factor. 
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